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CHAPTER I 
INTRODUCTION 


The diffraction by a, right-angled wedge is an important 
geometry in the analysis of many radiating systems. Although the 
diffraction of electromagnetic waves by a perfectly-conducting wedge 
has been studied extensively, few results exist for the case where 
one surface of the right-angled wedge is covered by a dielectric 
slab, an absorber coating, slightly rough or is corrugated. Such 
wedge configurations are of importance to problems of technology not 
only because they occur in practical structures but also because of 
design considerations where there is a need for the modification or 
control of edge diffraction. It is extremely difficult to find an 
exact solution to the diffraction by a wedge with coated or corrugated 
surfaces; however the problem becomes tractable if the boundary con- 
dition at this surface can be approximated by a surface impedance. 


A one-side dielectric loaded right-angled wedge is depicted in 
Fig. 1. If the dielectric slab is thin enough and the radiation from 
the end face of the slab at the edge of the wedge can be neglected, 
the dielectric-covered ground plane can be approximated as a reactive 
surface. This approximate electromagnetic boundary condition depends 
on the polarization of incident wave. For the TE case, the surface 
admittance, y $ , is given as 


*s 


(Sd) cot(Sd) 
- J z 0 W~^ 


9 


where £d is obtained from the solution of the transcendental equation 



OBSERVATION 

POINT 



Fig. 1. The 2-D conducting right-angled wedge with one side 
covered by a dielectric slab and a line source in 
the cylindrical coordinate system. 


(Cd)cot(5d) + j (e r -l)(kd) 2 -(Cd) 2 = 0; 

\ 

whereas for the TM case, the surface impedance 

z s = j -(ray x tan(5d) - 

i <* 

where (£d) is determined from another transcendental equation 
Ud)tan(gd) - \ (e f -l )(kd) 2 -(5d) 2 = 0. 


k and z Q are the wave number and characteristic impedance of free 
space, and d and e r are the thickness and the relative permittivity 
of dielectric cover. Admittedly the surface impedance is an approxi- 
mate boundary condition of limited validity, but its use may be ade- 
quate to approximate the important characteristics of the scattered 
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field; moreover, in principle it can be extended and improved[ll] . 

Other physical configurations which can be approximated by this 
impedance model are a conducting surface coated with an absorbing 
material, a special case of a rough surface[23] and a corrugated 
perfectly conducting surface with rectangular grooves, providing the 
groove spacing is small compared to the wavelength. 

The problem considered in this report is the diffraction of 
plane and cylindrical electromagnetic waves normally incident on the 
edge of a right-angled wedge; where one surface of the wedge is repre- 
sented by a uniform finite surface impedance and the other is perfectly 
conducting or has a surface impedance of zero. The diffraction of an 
obliquely incident scalar spherical wave is also treated. The solution 
is formulated in terms of two scalar fields which satisfy a uniform 
impedance boundary condition at one surface of the wedge and a Neumann 
or Dirichlet condition at the other. 

The difficulty of this mixed boundary condition can be overcome 
by a simple functional transformation. The transformation is made in 
such a way that the mixed boundary condition becomes a simple homo- 
geneous one. The Green's function resulting from the transformation 
can thus be found in the usual way. The final solution is obtained by 
an inverse transformation. This idea was first given by Lewy[15] and 
Stoker[24] who studied problems in water wave theory. The same idea 
has been applied by Karp and Karal[7,9], Chu[l], Karalet al[8] and 
Chu et al[2] for solving electromagnetic diffraction problems for a 
right-angled wedge; one of its surfaces sustains a surface wave, while 
the other is a perfectly conducting surface. The wedge is illuminated 
by a normally incident plane wave, a magnetic line source on the 
vertex or. an incident surface wave. In all the above cases, the 
solution is obtained by requiring that after the inverse transfor- 
mation, the derivative of the total field and the field itself must be 
continuous across the positive y-axis: x=0, y>0 (see Fig. 2). The 
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Fig. 2. The surface impedance model of a 2-D 
right-angled wedge with one wedge 
surface covered by a dielectric slab. 


solution is given in the far zone. In the case of a normally incident 
plane wave illumination, the diffracted field becomes singular at the 
reflection and shadow boundaries. Karp[10] also obtained a two 
dimensional Green's function for a right-angled wedge under an im- 
pedance boundary condition which does not support surface waves. One 
of the constants in his solution is, in general, difficult to compute. 
Using a different method we have obtained a simple expression for this 
constant. 
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In this report the transformation technique is extended to solve 
the diffraction of spherical waves. Also the relationship of the 
scalar solutions to the EM problem are considered. The transformed 
solutions consist of some special Green's functions which satisfy 
Neumann or Dirichlet boundary conditions at the surfaces of the wedge. 
Integral representations are obtained from the eigenfunction expansions 
of these Green's functions. The integrals are then evaluated asymp- 
totically via the modified Pauli -Cl emmow method of steepest descent[20] 
to obtain a far-zone approximation. The solutions are given in terms 
of a geometrical optics field and a diffracted field. The geometrical 
optics field consists of the incident field and the reflected field 
which comes from pole contributions to the integral. Saddle point 
contributions from the integral yield the diffracted field. A dif- 
fraction coefficient is obtained which is also valid in the transition 
region of the shadow and reflection boundaries. The geometrical optics 
field is discontinuous across the shadow and reflection boundaries; 
however the total field, which is the sum of the geometrical optics 
field and the diffracted field, is continuous. 

The scalar diffraction coefficients are derived in Chapter II. 

The wedge illuminated by a line source is treated first. By letting 
the line source recede to infinity, the diffraction coefficients for 
plane wave illumination are obtained. The problem for a point source 
illumination of the wedge can be related to that of a line source 
illumination by a Fourier transformation on the z (edge-directed)- 
coordinate. The scalar diffraction coefficients are composed of 
trigonometric functions and Fresnel integrels which are easy to com- 
pute. 


In Chapter III, the diffraction of normally incident electro- 
magnetic waves is treated. A ray-fixed coordinate system is introduced 
to represent the diffracted electromagnetic field. When the components 
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of the, incident and diffracted fields are expressed in this coordinate 
system, and the only leading term in the asymptotic solution is re- 
tained, the resulting expression for the dyadic diffraction coefficient 
reduces to a sum of two dyads; in the matrix notation, this diffraction' 
coefficient is a 2x2 diagonal matrix. 

It is shown in Chapter IV that the solution derived in the 
preceding chapters can be reduced to that of a perfectly-conducting 
wedge[20] and that it can be applied to a right-angled wedge with an 
impedance boundary condition illuminated by a magnetic line source 
at its vertex[9]. The radiation from a magnetic line source located 
close to the reactive surface of a right-angled wedge and the radiation 
from slots in truncated, dielectric-covered surfaces are considered. 

In both cases, the calculated and measured patterns are in good 
agreement. Patterns are given for different impedance-loaded wedges 
illuminated by a magnetic (or electric) current line source at dif- 
ferent locations, and some observations are made concerning the v 

behavior of the scattered field. nu. 

We conclude this chapter by briefly describing the related work 
of Felson[3], Malyughinetz[17,18] and Mohsen and Hamid[19]. Fe-lson 
has solved the scalar wedge diffraction problem for an arbitrary-, 
angled wedge when the surface impedance is proportional to the radial 
variable p. Because of this special boundary condition, the problem 
is separable and a Green's function can be obtained by standard 
methods. The uniform impedance boundary condition in our problem 
makes it impossible to apply the separation of variables method 
directly in this way. 

s • . ' ‘ ' i 

Malyughinetz considered the scalar diffraction by an arbitrary- 
angled wedge with different impedance boundary conditions on it s.^ two. 
surfaces; he restricted his solution to plane wave illumination. The 
solution is obtained by a method analogous to the one used by 


6 



Peters[21] and Senior[22]. With this method the differential equation 
and boundary conditions are expressed as a difference equation for the 
determination of a regular function whose real part represents the 
velocity potential . The exact solution thus obtained can be applied 
to our problem for the plane wave illumination if one of the surface 
impedances vanishes; however the integrals in his solution are diffi- 
cult to compute because they are improper with meromorphic functions as 
their integrands. At shadow and reflection boundaries, his solution 
appears to be- singular. 

Mohsen and Hamid have made a different approximation in their 
approach to the diffraction of electromagnetic waves by a perfectly- 
conducting arbitrary-angled wedge covered by a thin dielectric slab 
on one surface. They treat the diffraction of an E-pol ari zed plane 
wave normally incident on the edge. It is assumed that the dielectric 
slab is thin with the relative permittivity not much larger than 
unity and that the slab is in the illuminated region. By imposing an 
approximate boundary condition, E s 'v TE 1 , which is valid only in the 
far zone, at the dielectric covered surface, where r is the plane 
wave reflection coefficient for a dielectric slab on an infinite 
ground plane, the solution is obtained in an integral form, which is 
then evaluated asymptotically in terms of Fresnel integrals. 


CHAPTER II 

SCALAR WAVE DIFFRACTION 


In this chapter the diffraction by a right-angled wedge is 
formulated mathematically as a boundary value problem. One side of 
the wedge surface satisfies a uniform impedance boundary condition; 
the other a Dirichlet or Neumann boundary condition. Plane wave, 
cylindrical wave and spherical wave illuminations of the edge are 
considered separately. 
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A. Cylindrical Wave Illumination 


Consider a z-directed, uniform line source of unit strength 
radiating in the presence of a right-angled wedge with its edge also 
oriented in the z-direction (see Fig. 2). The total field U a , which 
consists of the incident and scattered components, satisfies D the scalar 
wave equation 

(1) (vj + k 2 ) U a = - fi(x-x') s(y-y') 

z b 

with boundary conditions 



8U a 


(2) 

9 V + *a^a 0 , <j> - 0 

" b b 


(3a) 

u a = 0 ; ♦•r ■ 


or 



(3b) 

3U b _ o . . _ 3ir 

H " ° » ♦ “ 2 


(4) 

Lim Sp (U p>a + j k U a ) -*0 , 0 < 4> 

<T- 


tt > 0 ° b b 

and 



(5) 

ll a is finite except at the point (x',y' 



b 

Here v t is the two dimensional Laplician operator, 6(x-x') is the Dirac 
Delta function, k is the wave number of the linear, homogeneous, iso- 
tropic medium surrounding the wedge. Condition (2) is a surface 
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impedance boundary condition; the subscripts a and b are used to 
denote the different surface impedances required in the EM problem 
to follow. In an acoustic problem, the surface impedance expresses a 
ratio of pressure to velocity; in an electromagnetic problem, it is 
a ratio of the tangential component of the electric field to the 
tangential component of the magnetic field. An impedance boundary 
surface may arise from an absorbing coating or thin dielectric slab 
on the surface or corrugations in the surface. Condition (3a) repre- 
sents a soft boundary for the acoustic case whereas (3b) represents a 
hard boundary for the acoustic case; mathematically these are the 
homogeneous Dirichlet and Neumann boundary conditions, respectively. 
Condition (4) represents the Sommerfeld radiation condition and (5) 
represents the Meixner edge condition. A time dependence, e Ja)t , 
assumed and suppressed. 

1 

Let us make a functional transformation 

(6) V = (|y + A) U. 

For simplicity, let U represent either U a or U^. By using the fact 

2 2 3 3 

that the operators (v. + k ) and (t- + x) commute and <s(y-y‘) = - 
3 toy oy 

<5 (y-y ' ) , V must be such that 
°y 

(7) (v* + k 2 ) V = (g^r - X) 6 (x-x 1 ) S (y-y ' ) 

and also satisfies the boundary conditions 

(8a). V = 0 ; $ = 0 

and 
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(8b) 



and the Sommerfeld radiation condition. It is worthwhile to mention 
that this transformation technique cannot generally be applied to an 
arbitrary-angled wedge since the boundary condition of V at the other 
wedge surface does not become simple except for some special angular 
space which, in this case, is a right-angled wedge. The function V 
need not be finite at the origin because it involves a differentiation 
of U. Let us consider a Green's function G which satisfies the 
conditions 


(9) (v* + k 2 )G = - fi(x-x') fi(y-y') , 0 < <j> < p- 

(10) G = 0 , <f> = 0 

and conditions (3), (4) and (5) for the wedge. It can be shown that 

00 

(11a) 6 = - i l e m J 2m (k^) H^)( kp> ) [ C0S ^ (*-♦') - 

m=0 -3 -3 

- cos— j (4> + <t>')]» for (3a) 

and 

(lib) G = - 1 l J 2mtl (k P< ) H<^(k P> ) [cos^SfU-*') - 

m=0 — ~T~ 

- cos-^y- (4>+4»’ )], for (3b), 

where • 
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o IK <f>' 1 j - 
0 1 P , p' < - , 

and 

e m = 1 , m = 0 

m 

= 2, m ^ 0. 

( 2 ) 

J (kp) and '(kp) represent the cylindrical Bessel function 
vm vm 

of the first kind and the cylindrical Hankel function of the second 
kind, respectively, and 


e< 



when p < p ' 

p > P ' 


p> = 



when p > p ' 
p < p ' 


Applying the operator - (gyr - A) to G, we obtain a unique 

function \p = - r «- x)G, which satisfies not only the conditions 

°y 

imposed in U, but also condition (5). This function is a particular 
solution of U. The complete solution to U is the sum of the particular 
solution and the homogeneous solution. By using separation of variables 
and applying the radiation condition, it can be shown that the homo- 
geneous solution is of the form 


( 12 ) 


v 

L 

m=0 






m 


cos v 


8 m sinv m *]. 


This homogeneous solution is singular. Since U must be finite, 
V cannot be too singular. Hence, of all the terms in Eq. (12), only 
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those terms for which v m < 1 are retained in the construction of V. 

Note a m = 0 and v m = if Eq. (3a) holds, while S m = 0 and v m = 
if Eq. (3b) holds. Thus V contains some parameters which are to be 
determined. 

We shall now introduce another approach to obtain the homogeneous 
solution for V[10]. Consider the Green's function W, satisfying the 
conditions (1), (3), (4) and (5), and the boundary condition 

(13) fy = 0 , * = 0. 

It can be shown that W can be given by a convergent eigenfunction 
expansion as 


(Ha) VI = - i l J 2m+ , (k P< ) (k P> ) [cos^ (♦-*•) ♦ 

m=0 3“ ”T" 

+ cos^y- (<f>+V ) for (3a). 

and 

(,4b) W - - $ Z «, J 2m (k P< ) «£>(*,,) [cos% (♦-♦• ) + 

m=0 ~3 ~3 

+ cos— (<J>+4> 1 ) 3 for (3b). 


From Eq. (1) we can prove that 


(15) 


<\ 2 + 


* 2 nff = + ly) ^x-x’) - 0. 


ft w 

Hence the function H e — - + — is a solution of the homogeneous wave 

oy oy 

equation. Furthermore, this function, H, obeys conditions (3a) (or 
(3b)) and (4), vanishes at 4> = 0, satisfies the radiation condition and 
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behaves at the origin like the non-vanishing terms in the series 
Eq. (12) which was used to represent the singular solution to be added 
to ill. Therefore, we obtain a complete solution for V as 


06) V- - (j|r- 1) 6+C(| 
where C is a constant to be determined. 


Once V is found, we can obtain the solution for the problem 
involving U since Eq. (6) can be integrated. A particular solution 
of Eq. (6) is 


(17a) 



V(x,?)d£ , 


Re X > 0 


(17b) 



Re x < 0. 


It should be noted that the above integrals are convergent and 
hence exist. In the far zone, we observe that U, V, G and W have the 
asymptotic formulas 


(18a) 

(18b) 

(18c) 

and 

(18d) 


G 


,-jkp 

,-jkp 


g(<t>;x' ,y') 


w -> — — - w($;x' ,y' ) 

Vp 

p-jkp 

V — v (<j> ; x ' ,y‘ ) 

VD 


U 


; -jkp 

IT 


U(4>;x' ,y ' ) , 
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Here G, V and W are known functions and U is to be found. Since 

the far-zone field is of interest, only those terms which are of 
*1 0 
0(~r) are retained, the operator — can be approximated then by 

-j k sin<j>; thus we obtain 


(19) (-j k sinet> + a)U - V, 

which can also be written as 


( 20 ) 


- jksin<j>W 
< C -H - j ks i n ' <|> 


+ XG - jksin<j>W 
x - j ks i n<|> 


In Appendix I, G and W are transformed into the integral repre- 
sentations and then computed asymptotically for the large parameter 
pj~r . From the asymptotic solution, we can identify the ray-optical 
behavior of the field. The incident field is 


U 1 = - i 




4 ^ Trk(p 2 +p ,2 -2pp ' cos (<t>-4> ' ))^ 2 


• » . 2 , 2 
-jk p +p - 


2pp 'cos (4>-4> 1 ) 


x e 


H (tt— | <#>-4> 1 1 ), 


where H(t) is the Heaviside unit function defined as 

(21) , H(t) {? ; JJ . 

From this we obtain 

(22) C = 1. 

Thus Eq. (20) becomes 


(23) 


II * - .iksin<j>W 

u x - jksin«j) * 
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The reflected field from the surface at <j> = 0 is given by 
( f N 


( 24 ) 


IT. = • R x 


-jk, 


V 


i _2J 

4 1 7Tk(p 2 +p ,2 -2pp'cos(<j)+<j)' )) 1/2 


p 2 +p ,2 -2pp 'cos (<{>+tf> 1 ) 


H (ir- | 1 I ) , 


where R is the reflection coefficient for an impedance surface, which 
is given as 


(25) 


n = jksincj>+A 
jksi n<j>-x 


3tt . 


The reflected field due to the perfectly-conducting surface <j> = -y is 


(26) 


■ (1 


ii 


4 iTrk(p 2 +p' -2pp 'cos ( ir -((()+<{) ' )) ^ 2 


-j k J p 2 +p ' 2 -2 pp ' cos (it- ) . 
x e H (tt— ) 3tt- (<(>+4> ' ) | ) 


where "-1" and "1" cire the reflection coefficients corresponding to 
the boundary conditions (3a) and (3b), respectively. The incident anc 
reflected field components are found from the pole contributions en- 
closed in the appropriate closed path which is given in Appendix I. 

The saddle point contribution yields the diffracted ray field. There- 
fore, in accordance with the geometrical theory of diffraction, U d can 
be written as 

(27) ut ^ U 1 D a 

b b /p 


in which 
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D is 

VS 

as 

(28) 

where 

(29a) 


(29b) 


(29c) 


(29d) 

and 

(30) 


u 1 ^ _ i Ji_ e -3^p ' 

U 4 irkp 1 e 

the diffraction coefficient for the boundary condition (3a) and 
the diffraction coefficient for condition (3b), which are given 


>a D s - Jfesin*. D q1 
a X - jksin<j> 

d 


“3 4 / + \ 

D , = - 1 2 sin£- [cos I- + cos |-]r + 

9 1 ZJZS i 3 3 3 J 

+ d*(s")F[Ka + (B")] + d‘(6‘)F[ica'(8")] + 
+ dg(6 + )F[«a + (B + )] + d-(s + )F[«-(B + )] 


d g< s > 


cos 


TT±S 

3 


sin 


1T±f3 

3 


D s = [d + (s')F[>a + (8‘)] + d'(6')F[Ka'(6')]] 
- jd + (8 + )F[i<a + (8 + )] + d“(6 + )F[Ka"(8 + )]] 


d ± (B) = - 




* 


. x b P q2 - J ksin * D h 
u b x b - jksin<j> 
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where 


(31a) 


4 r - + 

D 9 = " I 2 sin T [cosl cost-] + 

9 2 3/2^k t 33 3 

+ dg(e')F[Ka + (e")] + d'(8')F[ica'(0]}- 
- dg(e + )F[<a + (e + )3 + d'(s + )F[«V)] 


(31b) 


and 


(32) 


D h = )F[<a + (3 )] + d (B )F[<a (3 )]j + 

+ [d + (B + )F[Ka + ($ + )] + d'(B + )F[Ka‘( 3 + )]] . 



= ♦ ± <J>' , < 


kpp ' 

p+p 1 


F[<a (3)] = 2j Ka _ (B) 


e j*a + ( 3 ) 



dt. 


where the positive branch of the square root is taken, and 
• a ± (B) = 1 + cos(-S+3 N**). 

4 * 

The value of PT is determined by the integer which most closely 
satisfies the equation: 

3 N"tt - B = - ir. 


Outside the transition regions of shadow and reflection 

boundaries, (tea - ( b ) » 1. F[<a _ (3)] in Eq. (32) tends to 1. D$ Imd 

D a ] can be simplified as 
92 
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(33) 




D s = 
h 


2e ’ sin^- 
3 v® 


2tt r— 


1 


1 


[cosj^ - cos^-( 4> -4> ' ) cosj^ - cos|-(<j>- 


and 


(34) 


Dgl 

g2 


4e ^ sinj 
3 /2dT 


cos (^-) 


cos 


-) 


TT 


* _ 2 * 

[^OSj 1 - COS^U^' ) COSy 1 - COS^(<j>+<t>'_)j 


2 2 

which become singular when ^-( 4 >± 4 > ' ) = y^ + 2Nir, i.e., at the shadow and 
reflection boundaries. 


Thus, the region surrounding the wedge can be divided into three 
regions by the reflection boundary <f> = it - <j>' and shadow boundary 
<P = it + <j>', such as shown in Fig. 3.* The total field is given by 


(35) 


U * S 


U 1 + U r + U d 
U 1 + U d 


, 0 < <j> < IT - <j> ' 


, TT + (j) ' < (j) < 


3tt 

T~ m 


When the line source is near the impedance surface and far from 
the edge of the wedge, the surface wave is excited if the impedance 
surface supports a surface wave. The strength of excitation of the 
surface wave is assumed to be the same as those due to a uniform line 
source located above an impedance surface of infinite extent. Let 
l) s be the excited surface field. It can be shown that 


*We have chosen the example where only the surface 4>=0 is illuminated; 
the discussion for the two other possible cases follows in an 
analogous way. 


18 



OBSERVATION 
POINT P 



Fig. 3. The reflection boundary and shadow boundary of a line 
source illuminating on a 2-D right-angled wedge. 

(36) U = - e" X(y+y ' ) e' J ’ 3x , x>0 

5 P 

r~2 2~ 

where e = ik +x is the wave number of the surface wave. This 
excited surface wave propagates along the impedance surface <j>=0 and 
does not radiate until the impedance surface is terminated by the 
edge of the wedge. The source is assumed to be sufficiently removed 
from the edge so that the infinite surface representation is valid. 
Chu[l ,2] studied the radiation field of the surface wave diffraction 
by a right-angled wedge in which TE and TM surface wave diffraction 
were considered. TE surface wave diffraction corresponds to our 
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A q and B-j are constants and are given by 
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I is the Laplace transformation of H^(k|c|) and 



where -1 < Rev < 1 . 

\ 

In summary, the solution is obtained asymptotically in the 
far zone. When the impedance surface does not support a surface wave, 
the total field is given by Eq. (35). If the impedance surface 
sustains a surface wave and the line source is located close to the 
impedance surface and far from the edge, the diffracted field due to 
the surface wave excited by the line source must be taken into account. 
Note that the strength of excitation of the surface wave is pro- 
portional to e -Ay . Thus the excited surface wave has an important 
effect only when y' is small (<!>' « j ) . The total radiation field is 
then given by 

U 1 + U r + U d + V* , 0 < <t> < tt - ♦' 

(42) U * < U 1 + U d + ,ir - <J>' < <j> < tt + <f>' 

,ir + 4> 1 < 4> < 

Note that there is also a surface wave, C 2 e" A ^ + ^ x , propagating 
away from the edge due to the termination of the edge of the wedge. 

The strength of this surface wave, can be determined by another 
method which will be mentioned in Chapter V. This surface wave field 
is significant only when y is small (<J>1?0) and it does not radiate 
unless the impedance surface is truncated. The total field is the 
sum of U given by Eq. (42) and C 2 e “ A ^ +jex which exists only in the 
region {x>0, y>0}. 
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, The diffracted field, U^, due to the excited surface wave is con- 
tinuous in the whole region. The diffracted field, U , yields a 
uniformly asymptotic representation for the total field so that it 
properly compensates for the discontinuity in the geometrical optics 
fields, the incident field U 1 and the reflected field U r , across 
the transition boundaries, thereby, yielding a total field which is 
continuous everywhere. 

B. Plane Wave Illumination 


The plane wave incidence case can be treated as a special case 
of the cylindrical wave illumination. By letting p'-* 00 in the line 
source case of Section A and factoring out the "line source factor" 





2j_ o-JkP 


TTKp 


> 


the result for the plane wave case can thus be obtained. The details 
of the derivation are given in the Appendix I. Note that the total 
field, Ug, due to a plane wave normally incident on the edge of the ' 
wedge satisfies the 2-dimensional, homogeneous scalar Helmholtz 
equation 

(43) (v 2 + k 2 ) (p,W ) = 0, 

1 b 


where the subscript a (or b) indicates that the total field satisfies 
boundary condition (3a) (or (3b)), in addition to (2), (4) and (5). 


For a plane wave of unit amplitude and with its phase reference 
located at the edge, the diffracted field is given by 

d e" jkp 

(44) u° * D a i— - 

b b /p 
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where D a is given in Eq. (28) or Eq. (30). In the diffraction coef- 
ficients, the large parameter ^ i = kp since p 

C. Spherical Wave Illumination 


Consider the case of a scalar point source illumination of the 
wedge. The point source generates scalar, spherical waves. Let the 
point source be located at s 1 (p 1 ,<f>' ,z 1 ) as shown in Fig. 4. For a 
source of unit strength, the total field U a (p ,<f>,z;p 1 ,<P ‘ ,z ' ) due to 
spherical waves incident on the wedge satisfies 

(45) (v 2 + k 2 )U a (?,r ) = - 6(x-x')6(y-y' )s(z-z'), 

b 

and conditions (2), (4) and (5) as before, II satisfies the additional 

a 

boundary condition (3a), whereas U. satisfies (3b). Note that 

2 2 3^ D 

v = v t + is the 3-dimensional Laplician operator. Because of 
the special geometry of the wedge, this 3-dimensional problem can be 
reduced to a 2-dimensional problem via a Fourier integral trans- 
formation. Let 


(46) 


U a (p ,<f> »h ;p 1 ,<j> ' ,z ' ) = U a (s ,s ' ) e' jhz dz 
b -» b 


(47) U a (s,s 1 ) = U a ( P , 4 > ,h;p ' ,<j>' ,z ' ) e jhz dh. 

b * b 

Equations (46) and (47) form a Fourier transform pair. Let 
us assume the existence of the transform pair; then the Fourier trans- 
form of Eq. (45) yields 

(48) (v 2 +k^) U a (p ,4>,h;p ' ,4> ' ,z ' ) = - 6 (x-x 1 )<S (y-y ' ) e'^ 2 ’ 
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SIDE VIEW 


Fig. 4. Geometry for a point source illumination 
of a right-angled wedge. 
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where 


A - 


The z-variation of the function U a is removed, and rewritting the 
above equation as 


(49) 


(v^) U a ( P ,<f»,h;s 1 ) e 
b 


• c • i J hz ' = . 5 (x-x 1 )5(y-y' ) 


^ _ ihz' 

allows one to identify U a (p ,4>,h;s' )e J as U a (p,p') which is a 
2-dimensional scalar line source Green's function of (1), with the 
exception that k now is replaced by in (1). Hence U a (p ,<J>,h;F' )e^ hz 

is replaced by (p",p~' ;k t ). The total field U a (i"' ,£"' ) is thus obtained 

6 % b _ ihz' 

by taking the inverse Fourier transform of U a (p ,<f>,h;s' )e J . It will 

be shown in Appendix I that 


(50a) 


(50b) 


Ul = 6 


-jk | s-s ' 


U r = R - 


4-rr I S-S 1 I 
-jk|s.-s'| 


4tt|s.-s ' 


and 

(50c) 


U» * u’ 

§ slnB o 


s 1 -jks 
s(s+s ' ) e 


In which 


U 1 * S 


-jks' 


4irS 


r > 


s.. is the image point of the source s' and 
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jksincfr+A i 

I jksi n<j> -X , 

which depends on the boundary condition at the surface of reflection. 

Dg is again given by Eq. (28) or Eq. (30) with the exception that k is 
b i p l , 

replaced by sin S The ■ j -- ■ factor in Eq. (50c) is to be ex- 

o 

pected because of the conical spreading of the diffracted rays due 
to oblique incidence from the point source. This conical spreading 
is a consequence of the generalized Fermat's principle .[12] 

So far we have solved for the fields diffracted by a right- 
angled wedge due to different kinds of wave illumination^ we observe 
that the solution of the problem of a cylindrical wave illumination 
is essential for the cases of plane wave and spherical wave illumin- 
ation. The plane wave illumination is treated as a special case of 

cylindrical wave illumination when the source is removed to infinity 

, - • •• ^ 

whereas the spherical wave illumination is related to the* cylindrical 
case by applying the Fourier transform with respect to the z-variable. 

The asymptotic approximation in each case in the far zone yields 
a solution which can be interpreted in terms of the geometrical 
theory of diffraction. The pole contributions in the asymptotic 
approximation give the geometrical optics field which. consists of the . 
incident field and the reflected field, whereas the saddle point 
contribution gives the diffracted field. The geometrical optics field 
is discontinuous when one crosses a shadow or reflection boundary. 

This discontinuity in the geometrical optics field. is compensated by 
the diffracted field. Thus the total field, which is the sum .of 

• ' \r 

the geometrical optics field and the diffracted field, .is .continuous 
in the whole region. 
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It is interesting to note that the diffraction coefficient has 
the same form for the different wave illuminations. The only differ- 
ence for the different types of illumination is in the large parameter 
in the asymptotic approximation. If we let k = kL, then 


(51) 



for plane wave illumination 


L = ^ 


pp 1 
p+p' 


ss 1 
\ s+s 


2 

r sin 3 


o 


for cylindrical wave illumination 
for spherical wave illumination. 


In the next chapter it will be seen that the scalar diffraction 
coefficients D g and appear in the electromagnetic problem for 
plane or cylindrical waves normally incident on the edge. Furthermore 
it will be seen that when a ray-fixed coordinate system is introduced, 
the dyadic diffraction coefficient is merely the sum of two dyads 
or it can be expressed alternatively as a diagonalized 2x2 matrix. 
These properties were noted earlier in the asymptotic solution of the 

, • r i 

diffraction by a perfectly-conducting wedge[13,20]. 


CHAPTER III 

ELECTROMAGNETIC WAVE DIFFRACTION 


In Chapter II we treat the diffraction of a right-angled 
wedge illuminated by plane, cylindrical and spherical waves. The 
plane waves and the cylindrical waves associated with the uniform 
line source are normally incident on the edge of the wedge; the 
spherical wave, obliquely incident. In this chapter we show that the 
solutions of these scalar problems can be directly related to the 
corresponding electromagnetic problems provided that the electro- 
magnetic waves are normally incident on the wedge. 
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First, let us introduce the ray-fixed coordinate system des- 
cribed in Reference [13]. As mentioned earlier, the use of this 
coordinate system will enable us to write the dyadic diffraction 
coefficients as the sum of two dyads or in the form of a 2x2 diagonal 
matri x . 

Let us consider the field radiated from the source 0 and 
observed at P, as shown in Fig. 5. It is well known that the ray 



Fig. 5. Geometry for the ray coordinate system 
for electromagnetic wedge diffraction. 


incident on the edge at Q E in Fig. 5 gives rise to a diffracted field. 
To account for this Keller[12] has introduced a class of rays which 
includes an edge point such as Q E on their trajectory. Applying his 
generalization of Fermat's principle, the distance OQ^P along the ray 
path is an extremum, and the law of edge diffraction results. 
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Let the coordinates of the source 0 be (p 1 ,4> 1 ) and the obser- 
vation point P be (p ,<{>)» as shown in Fig. 5. The unit vector d is in 

A 

the direction of diffraction; the unit vector I is in the direction of 
incidence. The unit vector e is the unit vector parallel to z-axis 
(the edge of the wedge). Note that p = e x I, p^ = e x d, and 
e = z. The orthogonal unit vectors associated with coordinates 
(p 1 ,4> 1 ;p ,<j>) form a ray-fixed coordinate system. It will be shown 

A A 

later that when the incident field is decomposed into e- and p-directed 
components, and the diffracted field into e- and p^-directed com- 
ponents, the diffraction coefficients can be expressed as a 2x2 
diagonal matrix or as the sum of two dyads. 

In the following sections the diffracted field components will 
be determined first in the edge-fixed coordinate system, (p 1 , 4 > ' ;p ,<j> ) , 
which is the cylindrical coordinate system with z-axis aligned in the 
direction of the edge of the wedge, and then transformed to the ray- 
fixed coordinate system. 

A. Cylindrical Wave Illumination 

The diffraction of a normally incident cylindrical electro- 
magnetic wave by a right-angled wedge is considered in this section. 

The direction of incidence, which is normal to the cylindrical wave 

A A A A A 

fronts, is given by the unit vector I. Note that I • z = I '6=0. 
Cylindrical waves can be generated by a uniform electric current line 

A A 

source zl or uniform magnetic current line source zM located at 

(p 1 ,4> ' ) , where I and M are the current amplitudes of the uniform line 

sources. 

It can be shown that the axial incident electric field due to 

A A 

zl and the axial incident Magnetic field due to zM are 
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(52a) 


i-y 1 * c- iH^ckP-)] 


and 


(52b > h 1 *J 5 T M x [ - i H i 2 ) Ckp') 3 , 


respectively. The primes are used to denote the source coordinate 
system. 


First let us consider the cylindrical wave behavior due to the 
uniform electric current line source. The fields due to this line 
source are 

(53a) = 0, 

(53b > E z = -far 1 c ' { H o 2) (*!>')]• ' 

(53c) = 0, 

and 

(53d) ■ - ♦' ^ ^ 


( 2 ) 

For large kp ' , the Hankel function y (kp ■ ) can be replaced by 
its large-argument approximation. Thus 
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The equation kp ' = constant, describes a cylindrical surface. The 
direction of incidence 

(55) I = - p\ 

and 

(56a) I • E 1 = 0, 

(56b) I • H 1 = 0; 

hence in the far zone, the EM field is a cylindrical wave with its 
field vectors tangent to the wavefront, and perpendicular to its 
associated rays. The EM field of a uniform magnetic current line 
source can be shown to have the same properties. 

The total axial fields (E z and H z ) excited by uniform electric 
and uniform magnetic current line sources can be shown to be 

< 57a > E z=J^ IU a<^') 

and 

< 57b > H z = k MU b<^'>- 

respectively. The special function Ug (p”,pT' ) is the same as that of 
Eq. (1). The solution of Eq. (1) is given in Eq. (23). By substi- 
tuting the integral representations of Green's function G and W, 
shown in Appendix I, into Eq. (23), it becomes 
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where 


(59a) F a ^»+»+' ) " 67j(\-jksin*) f cot( ^T“) 


- cot (if-) - jksin<j> [sin (if) + sin (if- )] + 


costal) cos 2 (^S_) 

+ — + 1 — 

sin (if) sin( -j 8 ) 


(59b) 


F b (s -+*»' ) ■ ^ T x-jksinal i* I 51 "* 1 ?-) " s1n(^-) + 


4«/C+e 


cos 2 (if) cos 2 (if) 

+ — — - jksi n<#> x 

sin(if ) sin(if ) _ 


x [cot (if-) + cot (if-)] r . 


Note that b ± = <j> ± <j>' , f(s) = jO+cos£) and k = 

p+p 

The diffracted ray contribution to Eq. (58) is known from the 
result of Eq. (A-36) in Appendix I. Thus the axial diffracted fields 
can be written as 


(60a) 


d i e ~^ p 

z z a r* 


(60b) 


H z * H z D b 

z z b x 
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where E z and H z are given by Eqs. (52a) and (52b) with the Hankel 
function replaced by its large argument approximation, and D g and 
are given by Eqs. (28) and (30), respectively. 


For a problem which possesses cylindrical uniformity in 
structure and fields, the transverse fields E" t and can be de- 
termined from the longitudinal fields E z and H z using the following 
relations 


(61a) 


V. E 


E t = 


- jk z — | + jwu 


z X v t H z 


. , c ,2 ,2 

k -k z k -k z 


and 


(61b) 


V 


- jk 


v t H z 


. n 't E 2 
- Jwe o 9 

r-k : 


2 2 2 

where k = k t + k z . In the case of a normally incident EM wave, 
k^ = k and k z = 0. 


To evaluate E^ and H^, it is first assumed that we can inter- 
change the order of integration and differentiation when z x v 
operates on the integral, Eq. (58). Since only those terms of 0(— ) 

/p 

are retained in the asymptotic evaluation of the integral, the 
approximation ^ p is made. Thus 


(62a) 


V. E_ * p El e - jk(p+p 


t z 


[ -jk F a (€,♦,♦’) x 
a-L' 3 


1 jirk(p Z +p |Z -2pp 'cosC) 1 ^ 2 


e Kf(5) d; 


and 
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(62b) V H 'v p e" jk ( p+pl ^ 


-jk F. (c »4> ,4> ‘ ) x 
L-L* . D 


J jrrk(p 2 +p |2 -2pp ’COS?)^ 


e Kf(5) ds 


The diffracted field contribution to E z and H z is then 
obtained from the asymptotic approximation of the above integrals via 
the modified method of steepest descent with the result that 


(63) 


r * t *f 

7 t H z d 


* p (-jk) D a 


-jkp 




b So 


A A 


where D a in the above equation are identical to D a in Eqs. (60a) and 
(60b). Substituting Eq. (63) into Eqs. (61a) and b (61b), we obtain 
the transverse fields 

(64a) E 
and 

(64b) 


-jkp 


t ■ (E • p) p d u b -7= 

VO 


H t = (H 1 • p) p d D 


£jkP 

a 


—i * i 

In deriving the above equations, we have used H • p = - y z E' 

rj * ~ -j 0 ^ 

for a uniform electric current line source and E • p = z z H' for a 

o z 

uniform magnetic current line source. 


Thus the diffracted field for a uniform cylindrical electro- 
magnetic wave normally incident on the wedge is 

-i j ~ ~ 0 -jkp 

(65a) E^ = (E 1 ■ z)z 0 — 

and 
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(65b) 


tf 1 = (H 1 • p) p d D a 


A A 


,-jkp 


whereas the diffracted field for a uniform cylindrical electromagnetic 
wave normally incident on the wedge is 

-jkp 


(66a) 

and 

(66b) 


r 1 = (e 1 • p) S d D b 


■ g 

-jkp 


= (H 1 • z) z D b 

VO 


Let us define two dyads D E and such that 


(67a) 

and 


D e = e e D a + p p d D b 


(67b) 


d h = e e d k + P Ph D a > 


d a 1 


We can rewrite Eqs. (65) and (66) in the form of the geometrical theory 
of diffraction: 


(68a) 

and 


rd „ ri 


r^r(Q E ) • d e 


(68b) tid u i /n \ n e 


n 3 * H'(Q e ) • D h 


,-jkp 

"^T 

jkp 




; Alternatively, expressing the diffraction coefficients in matrix 
notation, we can write: 
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(69a) 


"c d 
^ H 


LPX. 


D a 0 


i_0 D b_, 


’ e !i(q e )" 

hlM 


,-jkP 

IT 


and 


(69b) 


,d 


'V 


D b 0 


•2 °a- 


h[(q e ) 

Lh;,(q e )J 


j-jkp 

IT 


where 


I r d - fd 
e, Ejl = E 


P d » Ej,(Q E ) = E 1 (Q e ) • e, and 
El(Q E ) = E i (Q E ) • p; Hj = • e, H*f, = H* • h]_(Q £ ) = H 1 (Q £ ) • e, 

• • A 

and H,|(Q e ) = H 1 (Q £ ) * p. Thus the matrix in the above equations is 
diagonal because the components of the field components are expressed 
in the proper coordinate system for the problem and only the terms 
of 0 (~t) are retained. 


B. 


Plane Wave Illumination 


For an electromagnetic plane wave normally incident on the edge 
of the wedge, the total electric field is given by 

(70a) E = z E 1 ’ U p 

d 

A i 

for an incident electric plane wave field z E , whereas the total 
magnetic field is given by 

(70b) H = z H 1 U p 

A i n 

for an incident magnetic plane wave field z H , where Ug is the 
solution of Eq. (43) and is given as ^ 
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(71) 


u? = 


L-L* 


F a (€,*.*') e jkpcosC dS 


with F a , 4 > , 4 > 1 ) given by Eqs. (59a) and (59b). 
b 

The diffracted ray contribution can then be obtained by 
asymptotically evaluating the integrals, Eq. (71), via the modified 
method of steepest descent. Thus, 

(72a) E^ * E^ D 

/p 

and 

a 0 “jkp 

(72b) H^H’D b l— , 

yp 

where E^ and are evaluated at Q £ , the point of incidence on the 
edge, and D a are given in Eq. (28) and Eq. (30) with the exception that 
k = kp in tfie large argument parameter. 


The transverse components of the diffracted fields can be com- 
puted in the manner similar to that of the cylindrical wave case 
treated in the last section. Following the same procedure, we obtain 

(73a) E* -v E i (Q E ) • D £ 

/p 

and 

(73b) FT -v H 1 (Q e ) * D h — 

v/p 


where 0 £ and 0 H are exactly the same as those given in Eqs. (68a) 
and (68b), except that k is equal to kp in this case. 
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Summarizing the results obtained thus far, the diffracted 
electric field can be written as 


(74) E^s) - E i (Q E ) • D E (s,I) A (s ) e' jks , 

where E -1 (Q E ) is the incident electric field at the point of dif- 
fraction, A(s) is the spatial attenuation which describes how the 
field intensity varies along the diffracted ray. 


(75) 


A(s) = — for plane and cylindrical wave incidence 
✓s’ 


D e (s,I) is the dyadic diffraction coefficient, 


(76) 


D c = zz 0 + 


PP d D b > 


where 

(77a) 

and 

(77b) 

where 

(78) 

and 

(79) 


_ x a D s - jksinf P q1 
a A - jksi n<}> 

a 

x b p q2 - jksin» D h 
u b A^ - jksinij) 

D s = |jjd + (g~)F[Ka + (0 )] + d (3”)F[i<a"(0~)]J- * 

+ -|d + (3 + )F[Ka + (3 + )] + d“(3 + )F[Ka~(3 + )]j-J , 

r - + 

Dqi = - 2 sin ?■ (cos I- ± cos %-) + 

g2 L 3 3 3 

+ {dg(3‘)F[Ka + (6')] + d-(e")FC<a"(0-)]} * 

* ^dg(3 + )F[Ka + (S + )] + dg(S + )F[Ka"(3 + )]Jj , 
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(80) 


-i- 

J 4 


d s (g) = - s— cot(4§.) 


3v 2rrk 


(81) 


♦ COS 2 (If.) 

V e) = sin (if) • 


in which 8 = 3“ = <p±<p‘ and 


(82) 


F[<a“(8 )] = 2j 


' Ka"(8) 


J<a ± (8) 


e-j* 2 dt. 


i Ka ± (8) 


The parameters which appear in F(<a ± (8)) are 


(83) a (8 ) = 1 + cos (-8+3 N ± tt) 

in which fT is the positive or negative integer whih most closely 
satisfies the equations: 


(84) 


3TrN + - 8 = ir 
_3irN - 8 = - t r . 


The quantity x = kL is the large parameter in the asymptotic evalu- 
ation of the pertinent intergrals involved in the formulation of the 
dyadic diffraction coefficient. The quantity L is given by 

for plane wave illumination 
for cylindrical wave illumination. 

When the ray-fixed coordinate system is employed, the dyadic 
diffraction coefficient can be expressed as the sum of two dyads. 

One of the dyads involves D . the scalar diffraction coefficient for 

a 
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TE (to 2 ) casej and the other, D^, for TM (to z) case. In turn, D g 
and are composed of trigonometric functions which are given by 
d ( b ) and d g (B ) and F(Ka"(e)) which involves the Fresnel integral. 

The latter can be regarded as a transition function which plays the 
most important role in the transition region of the reflection and 
shadow boundaries. The geometrical optics field is discontinuous 
across the reflection and shadow boundaries. It is this transition 
function which compensates for the discontinuity of the geometrical 
optics field and makes the total field continuous. Outside the trans- 
ition region, <a _ (e) > 10, F is approximately equal to 1. Curves of 
the magnitude and phase of F as a function of kLa are shown in Fig. 6. 
The scalar diffraction coefficients in the dyadic diffraction coef- 
ficient are independent of the wavefront of the incident field outside 
the transition region, as they are in the case of the perfectly- 
conducting wedge. 


CHAPTER IV 

APPLICATIONS AND NUMERICAL RESULTS 

The purpose of this chapter is to show that the diffraction 
coefficients and the analysis described in the preceding chapters 
are consistent with results obtained elsewhere for some special 
cases of our problem. The radiation from a magnetic line source 
located close to the reactive surface of a right-angled wedge and 
the radiation from slots in truncated, dielectric-covered surfaces 
are considered. In both cases, the calculated and measured patterns 
are in good agreement. Also some calculated patterns are presented 
to show the effect of surface impedance on wedge diffraction. We 
begin with the diffraction of an electromagnetic wave by a perfectly- 
conducting wedge. 
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Transition function F(kLa). 





1 . The Perfectly-conducting Wedge 


A dyadic diffraction coefficient for the perfectly-conducting 
wedge has been derived by Pathak and Kouyoumjian[20]. Our impedance 
surface reduces to a perfectly-conducting surface when x -*■-» and 

a 

x b = 0. In Eq. (28), it is quite obvious that when xy”- 00 , Dy^Dg , 
which is their soft wedge diffraction coefficient for the right- 
angled wedge, and in Eq. (30), when X^ = 0» ^ = which is their 
diffraction coefficient for a hard right-angled wedge. Thus the form 
of the dyadic diffraction coefficient, see Eq. (76), is the same for 
the impedance loaded wedge and the perfectly-conducting wedge, and 
so it follows that our solution reduces to that of the perfectly- 
conducting wedge[20]. 

2. Vertex-excited Surface Waves on one 
Face of a Right-angled Wedge 

When one wedge surface supports surface waves. Re x>0. For 

the line source placed precisely at the edge of the wedge, x'=y'=0. 

£ 

The first term in the series becomes singular. Thus the 

factor C-l , which is zero, multiplied by a singular term becomes an 
undetermined constant. Let the constant be C-| and Eq. (16) becomes 

(86) V = j | H| 2) (kr)cose + C-, Hj 2) (kr)cos |, 

3 

where e = - <f>. 

The above equation was also derived in a different way by Karp 
and Karal[9] with a different constant coefficient in the first term; 
their coefficient was given as jk/4. It can be shown that when 
x-*0, C^O. Thus V should reduce to the Green's function for a magneti 
line source located at the tip of a perfectly-conducting right angle 
wedge, which is H-^ 2 ^ (kr ) cos 0. For this reason we believe that 
our coefficient for the first term is the correct one. 
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3. A Magnetic Current Line Source Exciting 
a Right-angled Wedge 

In this section the radiation from a magnetic current line 
source located close to the reactive surface and far from the edge of 
the wedge is considered. The geometry is depicted in Fig. 7. The 


MAGNETIC CURRENT 
LINE SOURCE 



Fig. 7. Geometry of a line source radiating in 
.the presence of a right-angled wedge., 

reactive surface is such that.it supports a surface wave. Since it 
is assumed that the magnetic current line source is quite a few 
wavelengths from the edge, i.e., the exciting source is located in 
the far zone of the scattering from the edge, as noted in Section A 
of Chapter II, the total radiation f iel d ; i s the superposition of the 
geometrical optics fields and the diffracted field due to the rays 
incident on the edge and that due to the excited surface wave. Thus 
Eq. (42) is used to compute the radiation pattern. 

Chu[l] has considered this problem by using .reciprocity in 
conjunction with an exact solution and the, approximations involved in 
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the superposition. It is interesting to note that the expressions 
of the total radiation field are exactly the same except that the 
diffracted field, U d , due to the rays incident on the edge was not 
yet determined and neglected in his study. Figure 8 and Fig. 9 show 



Fig. 8. The amplitude pattern of a magnetic current 
line source exciting a right-angled wedge. 
(e/k=l .05 , L=10A, kh=l .5625) 
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Fig. 9. The amplitude pattern of a magnetic current 
line source exciting a right-angled wedge. 
(B/k=l .033, L=8.34X, kh=l . 309) 


0 ° 


the calculated patterns by Eq. (42) and those calculated and measured 
by Chu[l], The improvement by including the contribution of is 
evident in both cases. The ripple shown by Chu's computation in Fig. 

9 disappears because the contribution of U 1 ^ becomes more important 
when the magnetic current line source is closer to the edge. The 
radiation field pattern of a right-angled wedge with a reactive surface 
(8 ( k = 1.05)) and that of a perfectly-conducting wedge are shown in 
Fig. 10. The wedge is illuminated by a magnetic current line source 
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Fig. 10. The radiation pattern of a magnetic current line source. 
(L=l OX , kh=l .5625) 



located at L = 1 Ox and kh = 1.5625. The increase of the radiation 
field in the shadow region of the wedge with a uniform impedance 
surface is due to the fact that the energy stored in the trapped 
surface wave is radiated by the termination of the edge of the wedge. 

4. The Radiation from Slots in Truncated 
Dielectric-covered Surface 

Recently, Hwang et al studied the radiation from slots 
in a truncated, dielectric-covered surface[6]; in particular they con- 
sidered a narrow slot in the configuration depicted in Fig. 11. 



Fig. 11. The rays emanating from a slot in a 
dielectric-covered ground plane. 


Employing the geometrical theory of diffraction, the far-zone field 
is the sum of a geometrical optics field and a diffracted field. The 
geometrical optics field is the direct radiation from the slot to the 
field point P; it vanishes in the shadow region below the plane con- 
taining the surface AB. The slot also excites surface waves which are 
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incident at the termination of the dielectric cover at A and B. For 
simplicity the thickness t of the dielectric cover is restricted so 
that only the dominant TM Q surface mode is incident. The surface 
waves produce rays singly-diffracted from A and B to the field point 
P and to the lower edges C and D. They also produce reflected 
surface waves at these terminations. The contributions from the 
singly-diffracted rays and all significant multiply-diffracted rays 
are summed to give the field at P. 

The canonical problem for diffraction from edges A and B is 
the surface wave diffraction from a right-angle wedge, one of whose 
faces is perfectly-conducting, while the other is an impedance surface. 
The canonical problem of the diffraction of a surface wave on this 
right-angle wedge [2] is used to find the diffraction and reflection 
coefficients, which are adequate for sufficiently thin dielectric 
covers. However, the radiation from the vertical end faces AA' 

(or BB 1 ) of the dielectric cover cannot be neglected in treating the 
thicker dielectric covers, i.e., the radiation from the field inside 
the slab can no longer be neglected at its termination. The effect 
can be taken into account by a Kirchhoff approximation using the 
surface wave field along the vertical end faces of the slab. The 
problem is then treated as the radiation from a magnetic line current 
distribution on the perfectly-conducting surface of a right-angle 
wedge whose other surface is an impedance boundary. This contributes 
a second term to the diffraction coefficient previously obtained. In 
Fig. 11, the ground plane length L is 12.12 inches, the height H is 
1.71 inches, and the width is 22 inches. A 0.5 inch by 0.062 inch 
slot is positioned 0.1 inch to the right of the center of the ground 
plane, introducing a slight asymmetry in the pattern. The calculated 
and measured patterns are shown in Fig. 12; the two patterns are seen 
to be in good agreement. 
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Fig. 12. Pattern of a slot in a dielectric-covered ground plane. 
Frequency = 8 GHz. t = 0.1875 inch, e r = 2.56. 


5. Calculated Patterns of Impedance 
Loaded Wedges 

We conclude this chapter by presenting some numerical results 
which show the effect of the diffracted field and surface impedance 
on edge diffraction. 
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It is shown in Fig. 13 that the diffracted field properly com- 
pensates for the discontinuity in the geometrical optics fields, the 
incident field and the reflected field, across the transition 
boundaries, thereby yielding a total field which is continuous every- 
where. 


Figures 14 through 16 show the radiation patterns of a magnetic 
current line source illuminating a right-angled wedge. In each case, 
a magnetic line source illuminates the wedge from a distance p' = IX 
(wavelength) and at different 4>' locations. The normalized surface 
impedance, x^/k, is chosen to be 0, 0.1 and 0.5. The case x b /k = 0 
corresponds to that of a perfectly-conducting wedge. It can be seen 
from these patterns that when the reactance surface is in the il- 
luminated region, the side lobe in the vicinity of the reactance 
surface, which is the largest when the reactance surface becomes 
perfectly-conducting, is decreased, shifted and blended to the main 
lobe as x^/k varies from 0 to 0.5. Although the field may change 
considerably in the illuminated region far away from the reactance 
surface, it is almost the same in the shadow region for different 
surface impedances (note that this is not the case where the magnetic 
current line source is located close to the reactive surface which 
supports a surface wave). When the reactance surface is in the shadow 
region, the field does not change much as the surface impedance varies; 
except in the vicinity of the reactance surface, the field tends to 
decrease more rapidly as x b /k increases. 

The radiation patterns of an electric current line source 
illuminating a right-angled wedge are shown in Figs. 17 through 19. 

The normalized surface admittance, X /k, is chosen to be -1.0, -10 and 

a 

-°°. The case X /k=-°° corresponds to that of a perfectly-conducting 
wedge. It is noted that the radiation pattern does not change much 
as the surface admittance varies. 
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Fig. 13. Radiation pattern of a magnetic line source 
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Fig. 14. Radiation pattern of a magnetic line source 
Normalized surface impedance. 
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Radiation pattern of an electric current line source 
Normalized surface admittance. 





56 


Normalized surface admittance. 
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Fig. 19. Radiation pattern of an electric current line source 
Normalized surface admittance. 



CHAPTER V 
DISCUSSIONS 


It has been shown in the last chapter that the solutions for 
the radiation patterns of a line source in the vicinity of a wedge 
with a reactive wall are quite good for the analytic models con- 
sidered. As mentioned earlier, the effect of a coated surface which 
is covered by a dielectric slab or absorber is approximated by a 
uniform surface impedance. A uniform surface impedance boundary 

condition is then prescribed by + Xu = 0, where ^ is a normal 

an an 

derivative of u at the uniform impedance surface. In an EM problem, 
the value of X is given by either 



for a normally incident TE wave or 


x 



for a normally incident TM wave, where y g and z g are the surface 

admittance and the surface impedance of the surface, respectively, k 

and z Q are the wave number and characteristic impedance of free space 

and y Q = . The reflection coefficient can be shown to be 

o 


(87) 


d _ jksin<i> + X 
jksin<£ - X * 


In the case of a perfectly absorbing surface, the reflected field is 
zero. If the surface is perfectly absorbing at <f> = j which is in 
the direction normal to the surface, 

(88) X = - jk. 
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For a dielectric-covered ground plane with the thickness "d", 
it is easy to show that the reflection coefficient is given by 

jksi n4> + A a 

p - b_ 

j* " jksi n<j> - A 

the subscript 'a' denotes the TE case, ' b 1 denotes the TM case, 
(90a) x • = - 5 cot?d 

a 

(90b) A, = — tan^d, 

b r 

with 

(91) ? = k je^ - cos 2 <j > 

The total field, U a , satisfies the boundary condition: 

b 

(92) (f7 + ^a)U a = 0 

y b b 


(89) 

where 

and 


at the dielectric slab surface. Thus A a is a function of the 
variable 4> . We assume that if the dielectric slab is not very thick, 
A a can be approximated as a constant which is obtained by solving 
tne root, ?d, of the transcendental equation 

(93a) (?d) cotcd + |(e r - 1) (kd ) 2 - (cd ) 2 = 0 

for TE case and 


(93b) 


(ed) tancd - e » (e - 1 ) 


(kd ) 2 - 


(td ) 2 = 


0 
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for TM case. The root of the above equations corresponds to the 
surface wave mode if the dielectric-covered surface sustains a surface 
wave. It has been found numerically that the diffracted field is not 
very sensitive to the variation of the approximate surface impedance, 
but as far as the reflected field is concerned, this is not the case. 
Figure 20 shows the radiation pattern of an infinite, dielectric- 
covered ground plane illuminated by an electric current line source. 

The total radiation field is the sum of the incident field and the 
reflected field. In Fig. 20 it is apparent that there is a great 
deal of discrepancy between the exact radiation field and those 
computed by using the surface impedance model. This analysis, of 
course, shows an error in the reflected geometric field; A proportion- 
ate error will exist in the reflected diffracted field, i.e., the 
diffracted field associated with the reflection boundary, if the 
edge lies in the range 50°«}> 1 <130 o for this particular geometry. Out- 
side this range, the reflected geometric field should be reasonably 
accurate. The incident geometric and the incident diffracted will 
be little changed by the presence of the reactive wall. Thus, we 
suggest that when the dielectric-covered surface is in this illumin- 
ated region, the exact reflection coefficient in Eq. (89) is used 
instead of the approximate reflection coefficient obtained from the 
surface impedance model for the geometrical optics field. This would 
lead to a discontinuity of the order of 10 dB. The value of X should 
be taken from Eqs. (90a) or (90b) at 4> = in order to make the total 

field continuous across the reflection boundary. When the dielectric 
slab is very thin, it can be shown that 

(94a ) * a * - J- . 

and 

(94b) X. * — — k 2 d. 

b e r 

Thus the uniform surface impedance is a good approximation for this 
case. 


60 




ELECTRIC CURRENT 
LINE SOURCE 



S * 0.4 X 
€ r *>2.59 


Fig. 20. Radiation pattern of an infinite dielectric-covered 
ground plane (d = 0.127A, e r = 2.59) illuminated by 
an electric current line source (S = 0.4x). 


61 



The total solution of a differential equation is the sum of 
particular solution and homogeneous solution. The homogeneous 
solution of Eq. (6) is e _;Vy f(x). When A>0, r this solution is bounded 
and therefore exists in the domain of interest, which is recognized 
as a surface wave. Unfortunately our method, fails to find the ampli- 
tude of this surface wave. 

The strength of the surface wave was obtained for TM plane - 
wave illumination and a magnetic current line source located on the 
tip of the wedge by Karp and Karal[7,9]. The method they used is that 
after the functional transformation, the homogeneous solution of the 
transformed differential equation is chosen to be 

C 1 H{ 2) (kp) cos | 

3 

instead of that given by H {see Eq. (15)). The function 

C 1 h| 2) (kp ) cos | 

3 

does not satisfy all the conditions of the problem, since, following 
the inverse transformation its x derivative is discontinuous at the 
positive y axis. Thus the function e~ X ^ is introduced in the 
total field in the region xMD, y^O to make it possible to satisfy 
the continuity condition involving the x derivative of U. C-j and 
C 2 are then obtained. For the case of a magnetic (or electric) 
current line source located far from the edge of the wedge, the 
strength of the surface wave can be obtained in the same way. How- 
ever, two cases should be considered separately: and <f>' > j. 

In the former case, the line current source excites a surface wave 
which propagates toward and away from the edge. The strength of 
excitation of the surface wave is assumed to be the same as those due 
to a uniform line current source located above the same impedance 
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surface of infinite extent. The justification of this assumption is 
that the line source is far from the edge, i.e., the exciting device 
is in the far zone of the scattering from the termination which is, 
in our problem, the edge of the wedge. The surface wave field is 
well established and propagates along the impedance surface before it 
reaches the edge. Thus in matching the continuous condition at the 
line (x=0, y>0}, this excited incident surface wave should be included 
in the total field. However, in the latter case, there is no such 
incident surface wave included. 


Several approximations have been made in deriving Eq. (90). 
According to the geometrical theory of diffraction, it has been 
assumed that the high-frequency diffracted field propagated along 
its ray path perpendicular to the diffracted electric and magnetic 
fields in the same way as the geometrical optics field. This approxi- 
mation does not introduce a large error if both the source point 
and the field point are far from the edge. Second, in deriving Eq. 
(90), it has been assumed that kL is large. However, based on the 
extensive numerical study of asymptotic solutions of this type pre- 
sented in References [5] and [20], it would appear that this approxi- 
mation generally introduces serious error only when kL < 1 . 


The preceding discussion has been restricted to the diffraction 
by wedges with straight edges; however the geometrical theory of 
diffraction can be used to treat the diffraction from curved edges[20]. 
The diffracted ray paths are determined by the generalized Fermat's 
principle for edge diffraction, and the conservation of power flow 
in the resulting astigmatic bundle of rays, see Fig. 21, leads to the 
general spatial attenuation factor 


(95) 


A(s) = JsTpT+sJ 
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DIFFRACTED 
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where the caustic distance p shown in Fig. 21 is given by 



for spherical wave illumination of the edge [14]. In the equation 
above 

P e is the radius of curvature of the edge 
n is the unit vector normal to the edge 

A A 

I ,s are unit vectors in the directions of incident 
and diffraction, respectively, 

A Ak 

e is the angle between I, e, the tangent to the 
0 edge at the point of diffraction, see Fig. 22. 

The expression for A(s) given in Eq. (95) will reduce to that 
given in Eq. (75) if p. is set equal to infinity for the straight 

cr 

edge; furthermore, for plane and cylindrical wave illumination of 
the straight edge, p c = 00 . 

The diffraction coefficient is assumed to be independent of p g 
to a first approximation; this is a reasonable assumption because it 
is independent of the curvature of the incident wavefront to this 
approximation, as has been demonstrated in Chapters II and III. 
Furthermore, the validity of this assumption has been confirmed in 
numerical applications of the geometrical theory of diffraction to 
structures with curved edges of a perfectly-conducting wedge[14,20] . 
Thus, in accordance with the postulates of the geometrical theory of 
diffraction, Eq. (74) becomes 

(97) f(,> « e’(Q e ) • 5 e <s,I)L^ e-j ks 

c 
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where Dg is identical to Eq. (76) with = j. Alternatively, if the 
incident and diffracted fields are decomposed into components parallel 
and perpendicular to the plane of incidence and diffraction, respec- 
tively, we can write Eq. (97) in terms of matrix notation as 



By introducing the ray-fixed coordinate system the polarization 
effects of high-frequency scattering can be greatly simplified, 
whether this involves the reflection from a smooth curved surface, 
the diffraction from an edge, or the diffraction from a smooth curved 
surface. Specifically, the polarization of the scattered field can 
be related to the polarization of the incident field by a 2x2 diagonal 
matrix. 

Several important steps used to solve this boundary valued 
problem are briefly described in the following to show the merit of 
the method employed. The mixed boundary condition of a wave equation 
is transformed to a simple homogeneous one by a functional trans- 
formation. In general, the transformed wave equation becomes compli- 
cated at the expense of the simplification of the boundary condition. 

g 

By using the fact that the operator, L = — + A, of the functional 

ay 

2 2 3 3 

transformation and (v + k ) commute and gyS(y-y') = - gyr<5 (y-y ' ) , 

the solution of the transformed wave equation can be related to some 

special Green's functions which satisfy the simple homogeneous 

boundary conditions. Although the technique is applied for the 

right-angled wedge in this report, Karp[10] pointed out that it could 

also be used to simplify the boundary value problem for some special 

angular spaces. Currently, this method is extended to treat an 

arbitrary-angled wedge with one of its surfaces satisfying a uniform 
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surface impedance boundary condition while the other, a perfectly- 
conducting surface, satisfying Dirichlet boundary condition. The far 
zone approximation at the perfectly-conducting surface is made after 
the functional transformation so that Dirichlet boundary condition 
remains unchanged. Thus the boundary conditions of the transformed 
wave equation still become homogeneous. The method can also be 
employed for a right-angled wedge with two uniform impedance surface 
walls. Recently, the diffraction of a two-part plane by a plane wave 
illumination is already pursued by the Wiener-Hopf technique; one 
part of its surface satisfying a uniform surface impedance boundary 
condition while,, the other being perfectly-conducting. The trans- 
formation technique fails for this case since the boundary conditions 
of the transformed wave equation can not be simplified by using the 
functional transformation. 

In summary. Maxwell's equations can be solved exactly for few 
problems with appropriate boundary conditions. Hence approximations 
have to.be made to obtain solutions to useful problems . High-frequency 
asymptotic approximation is used to solve the diffraction problem. 
Mathematically speaking, by making a transformation which simplifies 
the boundary condition and employing the modified Pauli-Clemmow 
method of steepest descent, the diffraction coefficient has been 
found for the right-angled wedge with one wedge face satisfying a 
uniform impedance boundary condition. The introduction of a ray- 
fixed coordinate system yields a compact, dyadic diffraction coef- 
ficient for EM waves illumination. This diffraction coefficient is 
also valid in the transition regions of the shadow and reflection 
boundaries for a variety of edge illuminations. The geometrical 
theory of diffraction is used to extend the solution to more general 
cases, such as the curved edge. The diffraction coefficient is 
composed of trigonometric functions and Fresnel integrals which are 
easy to compute. 
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APPENDIX I 

AN ASYMPTOTIC APPROXIMATION FOR SCALAR WEDGE DIFFRACTION 
VIA THE MODIFIED STEEPEST DESCENT METHOD 


A. Cylindrical Wave Illumination 


Let us consider . a scalar field U a which satisfies the wave 
equation b 

(A— 1 ) (v 2 + k 2 ) U a = - 5(x-x') fi(y-y') 

* b 

with the boundary conditions 


(A-2) 


(A— 3 ) 


(gy + x a) u a 
dy b b 


= 0 ; <p = 0. 


r 


< 


U a 

3U fc 

dip 


= 0 


= 0 


<j) = nu 


and which also satisfies the Sommerfeld radiation condition and the 
Meixner edge condition and n=1.5. In Eq. (20) the far-zone form of 
U a was shown to be 


(A-4) 


U a * 
b 


[(C-l)a|T+ A a ] G - jkCsin*W 

3y b 

-jksin<Ji + X a 
b 


where G and W are two known special Green's functions given in 
Eqs. (11a), (lib), (13a) and (13b), and C is a constant to be 
determined. The geometrical configuration is depicted in Fig. 2. 
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An integral representation for the product, J v (kp')H^(kp) 
is Reference [16] 


m m 


12) ^ f c -j” 4ct-k 2 (p 2 +p ,2 )t _1 ] k 2 , 

(A-5) J (kp ' )H^ Z; (kp ) = - -U- e 2 I (K_££_) 

V m V m TTJ v m ' t ' 

mm ' u m 


k 2 pp 1 


where C > 0, v > - 1 and |p 1 1 < |p | . I u ( K ) represents the 
m m z 

modified cylindrical Bessel function of the first kind. The integral 

k 2 nn • 

representation for I v (- ~ r— ) is given by Reference [16] 

m 


(A-6) 


V * 


1_ 

2tt 


2 

Y«+joo K BP cos g + jv £ 

e 1 m d5 

Y+j“ 


where - tt < y‘ <0 and ir < y < 2ir. Thus 


(A-7a) 


i (^> ■ 

v ' t 
m 


1_ f 

2n 


ifi cos? + j V 

e d? 


or 


(A-7b) 


i (*V-> - 
\ 4 


k 2 pp 1 


1 

2ir 


cosC 


JVE 

m 


-'L' 


dC 


The contours L and L' are indicated in Fig. 22. Let 6 = 

then Eqs. (lib) and (13a) can be written as 

(A-8) G, = g({T,p"';B’) - g(p,p';B + ) 


and 


(A-9) W = g(p,p';B") + g(p,p';B + ) 

a 


dt 

t * 
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Fig. 22. The L agd L' contours for the integral representation 
of I m ( k ) in ^e complex s-plane. 
n 
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g (p t p ) - - 8j 2 n ) 0 


1 ,c-j- ^lpV 2 lf'] 


k 2 pp’ cose 50 jp£pU + e) 

# T l e 2n & + 


+ e 


. C I .OVT 

Mr- cos£ * "3 2n 

> 1 I e 


.1 dt 


d ?r 


, : . . ' cotef^) - k 

(A-12a) 1 e 2J Zn 


- -3^6) , ( ea ) + T. 

(A-12b) I e 2o • 2n 

m=0 


(A-13a) 


I e 

m=0 




2 / C + 3 \ 

, cos tap/ . p+$\ 

1 feD — + sin tap) » 


(A-13b) 


I e 

m=0 




2 /E+&\ 

COS (oT-i . P+Sn 

1 Q — + sin tap) 

7TT Tr + ft v £** 


2 3 sin(|p-) 


and Reference [16] 


Z 

2t dt 
t 


- K o (jZ) ' 
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Equation (A-ll ) can be simplified as 


(A-15) g(p,p;g) 


o 2 . 

8ir nj 


L-L' 


dC 


cos 


[ 2r\ } 


s ’ n <§r> 


+ 


K 0 (jz(e)) 


where K_ (JZ) is the modified cylindrical Bessel function of the second 
o I s 2 " 

kind, of order zero and argument jZ, and Z(0 = k^p +p ' - 2pp ' cos£. 

We are interested in solutions for which Z(0 is sufficiently large 

which will be justified later. 

If |Z(c)| is large, K q ( jZ ( c ) ) can be replaced by its large 
argument approximation 


(A-16) K o (jZU)) * J 23 ZTTT e " jZU) ‘ 


In the asymptotic solution described later, it will be seen 
2 2 

that p +p' -2ppcos£ >> 0 in the neighborhood of the saddle points, 
therefore 


(A-17) 


-jkZ(c) _ e -jk|(p+p ' ) 2 -2pp '(1+cosS) 


'v e 


■jk(p+p')[l - — (1+coss)] 

(p+p 1 ) 


Thus 

(A-18) 


{ - - 2 i ttj 

L-L 1 ^ jirk(p +p‘ -2pp'cosC) ' 


F^S.g) e Kf(5) dC 


where 
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(A-19a) 


= 4irnj 


cos^( 


5+3 

2n 


) 




+ sin(|^-)| 


(A-19b). f(5) = j(l + cos?) 

and 


(A-19c) 


K 


- kpp 1 
p+p' * 


Therefore the integral in Eq. (A-18) is in the proper form 
to be evaluated by the method of steepest descent for a large 
parameter. The saddle points of f(s) occur at 




= 0 ; 


but only 5 S = ± it are considered because the steepest descent paths 

through 5 s = ± -n allow us to close the (L-L ' ) contour. It is clear 

that in the neighborhood of these saddle points the inequality in 
2 2 

p +p' »2pp'cos is satisfied, which justifies the approximation in 
Eq. (A-16). Figure 23 shows the locations of the steepest descent 
paths through the saddle points at 5 g = ± ir. Therefore 


(A-20) 


,-jk (p+p ' ) 


= - e 


d5 F 1 (c,3) i ? 2 

L-L' jirk(p +p -2pp 'cos?) 

jk (p+p ' ) 
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+ 1 d? F, (5 ,3) s 2 T7? 

SDP(tt) 'SDP(-tt) 1 1 jtH<(p +p -2pp 'cos?) 

. r fThe residues of the integrand enclosed 
7TjS lby L-L' , * * 


:f(5) 


+ 2wj*-iu.. , ■ • and SDP ( ±7T ) 

+ Branch cut contributions, if any. 
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Fig. 23. Steepest descent paths (SDP ± tr) and 
the complex K plane topology. 


The pole singularities occur at 

C p = - $ + 2nNir , N * 0, ± 1, ± 2, ••• . 

The residues corresponding to are evaluated only for those |c ] <jr . 
Let the poles contribution to g be denoted by g p , thus 


(A-21 ) 
g p (p~ ,p~' ;e) 



2 -2pp'cos5 ) 1/2 

1 COSNtt X 


x H[tt- 1 6+2nNu | ] 
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where 


r 


0 , t < 0 

(A-22) H(t ) = 1 , t = 0 

1 t > 1 




The saddle points contribution to g, which is denoted by g d , 
is derived via the modified steepest descent method[20], where the 
pole singularity close to the saddle point. is taken into consideration. 

rl ' . * 

Thus g a can be written as 


(A-23) 

where 

(A-24) 


g d (p ,p 1 ;e) * 


i -jk P ' 

4 lirkp 6 ' 


p -J k P 

d — 

9 G 


-i- 
_ 3 4 

d 9 (e) " " 2n/HIT 

. c ° s2 <gr> 

Sin(^) 


sin(gS-) + sin(^) 


F[ica + (s)3 + 




F[ica"(e)] }■ 


J 


with 


(A-25) F[Ka ± (e)] = 2j ^a*(B) e jKa " (e) 


..2 

e" Jt dt, 


W( 3 ) 


where the positive branch of the square root is taken, 

a ± (3) = [1 + cos{- 3 + 2nN ± Tr)], 
and 


= k PP 

p+p ' 
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The value of N is decided by the integer which most closely 
satisfies the equation 

2nN~ir - 0 = ± it. 

Substituting Eqs. (A-21 ) and (A-23) into Eqs. (A-8) and (A-9), 
we obtain 

(A-26a) g£(p,p') = g p (p,p';e~) - g p (p,p';e + ) 

(A-26b) GjJfp'J') = g d (F,p';8') - g d (p _ ,p"' ;e + ) 

and 

(A-27a) W p (^,p') = g p (p,p';8~) + g p (p,p';3 + ) 

(A-27b) w d (^') = g d (F,p' , ;e") + g d (p,F';e + ) 

where the superscript p denoted the contribution from pole singulari- 
ties; the superscript d, from the saddle points. 

Equations (11a) and (13b) have been evaluated in the same manner 
by Pathak and Kouyoumjian[20], in which G a and W. are given as 

a D 

(A-28a) G p (p* »p"' ) = I p (p,p‘;0") - I p (p‘.p , ;0 + ) 

(A -28b) G d (p” ,p*' ) = I d (p\p* ;B") - I d (F,r';6 + ) 

and 

(A-29a) w£(p,p') = I p (p,p';6 ) + I p (p,p';S + ) 

(A-29b) wjjGT.p’') = l d (p'»F , ;0") + l d (p’,p 1 ;B + ) 
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where 


(A-30) I p (^’;6) = 




j7rk(p 2 +p ,2 -2 pp 'cos£p)^ 


-jk(p 2 +p ,2 -2pp 'cose ) 1/f2 
x e p 


x H [tt— -3+2nNir | ] , 


and 

(A-31 ) 
with 


(A-32) 


I d (p,p';3) * 




e -jkp 


d 


e~J kp 



The superscripts p and d of I(p,p';3) denote the contributions from 
pole singularities and saddle points, respectively. 


Now, let us consider the pole singularity contribution for the 
case N = 0 and 3 ” = The equation | 1 1 <it , i.e., <}><tt+<j>', 

describes an illuminated region for a cylindrical wave illumination 
on a wedge (see Fig. 3). The spatial factor 




p 2 +p ,2 -2pp 'cos U-V ) 


\ p 2 +p' 2 -2pp 'cos{4>-<}>' ) 


indicates that this pole singularity contribution yields the incident 
field. Thus, substituting Eqs. (A-26a), (A-27a), (A-28a), and (A-29a) 
into Eq. (A-4) and recognizing that U is the incident field. 
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i 2.1 -jk/fcvT 

'4 irkp 1 e 

we obtain 

C = 1. 


Therefore Eq. (A-4) becomes 


(A-33) 


U a (p ,p ' ) * 
b 


x a G a - jksin<f>W a 

b b b 

X a - jksin<}> 

b 


The above inference can be reverified by considering the 
reflected field, which is obtained from the pole singularity corres- 
ponding to the case N = 0 and & + = , we obtain 


(A-34) 


jksin<f>+X a 

gP = b 

| jksm<|>-x a 


. i 

2 .i 

4 i 

frk(p^+p '^- 2 pp 1 cos ( <f>+ 4 > ‘ )) ^ 


x e -jk(p 2 +p , 2 - 2 pp 'cos ( 4 >- 4> 1 )) 1/2 , 


x 


which is the reflected field from the reactance surface, <j>=0. Whereas 
for the case N=1 , e= <h*4> ' , 


(A-35) uj?(p ,p ' ) = 


-1 
L1J 


, _ i 2 j 

[_ 4 Trk(p^+p' -2pp 1 COS (2nir-(<j)+<j) 1 ) ))^ 

x e -jk(p 2+ p 1 2 -2pp 'cos (2mT-(cj>+<j> ' )) )^ 


is the reflected field from the perfectly-conducting surface, <j>=mr, 
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The diffracted field component is then obtained from the 
saddle point contributions to the function G and W, which is given 
in the form of 


(A-36) 


where 


a — • -j kp 

Ug(pV ) * U 1 D a , 

b b Jp 


U 1 ^ 


i, 

4 Ttkp 


e 


-jkp 


I 


(A-37a) 

(A-37b) 

with 


D = 


D l = 


x a °s - j ks1l "* p ql 
A - jksin 4 > 

a 


\ °q2 - jKsin^D h 
A^ - jksin# 


D s = d((T) + d(e + ) 
h 


and 



■ d 9 (s‘) ± d g (e + ) 


B. Plane Wave Illumination 


The result for the plane wave case can be obtained by letting 

p in the line source case of the last section A so that •> kp 

p+p 

in the expression for k. However, for the sake of completeness, we 
shall treat this problem by starting from Eq. (A-5) in which p and p 1 
are interchanged. In this case, kp ' >> kp . Thus, instead of the 
approximation made in Eq. (A-16), we have 
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(A-38) 


K Q (jkZ(c)) * 




2tt e jkpcos?: 


which is given in Reference [16]. Therefore Eq. (A-18) becomes 


(A-39) 


g(p,p';6) ^ 





F 1 (5,8) e J * kpcosC d5 . 


Eliminating the "line source factor" appearing in the bracket and 
proceeding in a manner similar to that of part A for cylindrical 
wave illumination, we obtain 


(A -40a) 


g p (p ,p 1 ;3) 


jkpcos? 

= e p 


cosNu x H[ir- [ -$+2nNir | ] , 


and 


(A -40b) 


g d (p,p';3) 'v d g 



where d g is identical to Eq. (A-24) except k = kp in this case. 

The asymptotic form for and W a are also given by Eqs. 

(A-26) and (A-27), respectively. G a and for the plane wave illum- 

ination are derived in the same way by Pathak and Kouyoumjian[20], in 
which they give 

_ jkpcosC_ • •• 

(A-41a) I p (p ,p 1 ;3) = e p H[tt-| -0+2nN-rr | ], 

and 

h -jkp 

(A-41 b ) I a (p ,p 1 ;6) -v d 

yp 

where d is the same as those given in Eq. (A-32) except the large 
parameter, ^-r , replace by kp. G p and G d are still given by Eqs. 
(A-28) and (A-29). 
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In summary, plane wave illumination on a wedge is a special 
case of cylindrical wave illumination. The diffracted field com- 
ponent for plane wave illumination is thus still given as 

H . _• -jkp 

(A -42) U°(p,p') * U 1 (Q ) D, , 

b E 5 


where D a are given in Eqs. (A-37a) and (A-37b) with the large parameter 
. , b 

replaced by kp . The reflection coefficient R is given as 


(A -43) 


R 


1 

-1 

jksi ncf>+x 
jksin<j>-A 


which depends on the boundary condition at the surface of reflection. 


C. Spherical Wave Illumination 

Consider a scattering field U a (F,F' ) which satisfies the wave 
equation b 

(A -44) (v 2 +k 2 ) U a (?J') = - 6 (x-x 1 ) <5 (y-y 1 ) 6(z-z'), 

b 

and the boundary conditions (2), (4) and Meixner edge condition, 

U g satisfies Eq. (3a); U b> Eq. (3b). Fourier transforming Eq. (A-44) 
yields 

(A-45) (v 2 +k 2 ) Qf a (p.*,h;p , ,*',z') = - 6 (x-x') 5 (y-y') e' J * hz 

x 1 b 

where Ug and ftg are a Fourier transform pair given in Eqs. (34) and 
(35), and k^ = k^ - h^. Equation (A-45) can be rewritten as 
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(A-46) (v?+k?) Of a (p ,(j) ,h;p ‘ ,<1> ' ,z 1 ) e jhz ' = - 6(x-x') 6(y-y*) 

t x b 


which is identical to Eq. (A-l ) with the exception that k is re- 
placed by k t in Eq. (A-l). Thus 


(A -47 ) (|y + X)tt a e jhz ‘= [(C-l)^-+ X]£ a (p ,p ' ;k t ) + |y \^g(p ,P ' ;k t ) 

where 

= ?(p,p';B ,k t ) + ?(p ,p ' ;0 + ,k t ) 

and 


(A-48a) 


l 


(A-48b ) 

with 

(A-49a) 

and 


X 


W. 


a- 1 


= g(p,p';6 ,k t ) + g(p , p ' ;$ T , k t ) 


— . + 


g(p ,p',0,k t ) 


1 

8tt nj 1 


L-L 1 


c°s 2 (fr-) 


sin(^) 


C+B, + S 1 ' n( fr> 


K 0 (jZ({))d5 


(A-49b) TW;&,k.) = -J- f cot(f^-) K (jZ(c))d5 

1 8/nj ■’L-L' 0 


2 2 

Z(?) is given as k^p +p' -2pp'cos£ . 

From the inverse Fourier transform relation and the approxi- 

£ 

mation, r— ~ - jksin^ , made in the far zone, Eq. (A-47) becomes 
°y 
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(A-50) 


U a (s,s') * 
b 


[(C ‘ 1 W r + Ga " jkcsin^Wa 

_ b b_ 

X - jks i n<#> 


G a and W a can be written as 
b b 


(A-51a) 

and 


L W b-i 


= I (p ,z;p 1 ,z' ;$") + I (p ,z ;p ' ,z ' ;B + ) 




(A-51b) 
where 

(A -52a) I (p",z ;p"' ,z' ;3) = 


= g(p,z;p' ,z';0 + g(p ,z;p 1 ,z ' ;g + ) 


1_ 

2-rr 

and 

(A-52b) 

gtp’.Zip"' »z 1 ;b) 


L-L 1 8tt jn 


~ cot( lr' )K o (jz{5))e ’ jh(z, ’ z) d5 


dh. 


1_ 

00 

f 1 f 

2tt 

• 00 

'L-L' 87T^jn [ 


cos ^ 

*1n(&) 




+ sin (ffp) 


KjjZU))e" jh(z '~ z, d£; 


y 


dh. 


By utilizing 


• ^ n 

(A-53) K v (jZ) = - j f- e ' 2 H< 21 (z) 


and[4] 


83 



(A-54) 


dh ^{p^p-^pp-cosc} 1 / 2 ) 


p -jkp 2 +p' 2 -2pp'cos5+(z-z' ) 2 

j j= _ = - = 

\ p 2 +p' 2 -2pp'cosC+(z-z' ) 2 


Equation (A-52b) becomes 
(A-55) g(p*,z;p’ 1 ,z‘ , 0 ) = 


FtU.6) 
L-L* 1 


e -jk, p 2 +p ,2 -2pp 'cos5+(z-z' ) 2 
4-m p 2 +p |2 -2pp'cOSS+(z-z' ) 2 


where F-|U,e) is given in Eq. (A-19a). 


Note p=s sine Q and p'=s' sin 6 o (see Fig. 4). Hence 
| (Z-Z 1 ) 1 1/2 = (s+s 1 ) cose o , and 

— 2 

o , 99 2 ss'sin e 

(A-56) p +p' -2pp'cos5+(z-z ' r = (s+s 1 ) 2 1 5 ^- x 

(s+s 1 ) 2 


x (l+cose) 


Substituting the above equations in Eq. (A-55) yields 


(A— 57 ) g(p ,z; P 1 ,z' ;g) 


» 

dC F,U,0) 
L-L 1 1 


2ss'sin^g 

-jk (s+s * ) 1 o— - ( 1 +cosC) 

J (s+s 1 ) 2 

4 tt, (s+s 1 ) 2 - 2 ss‘sin 2 B 0 (l+cosc) 
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As before, the integral over L-L' is evaluated asymptotically for 

ksc 1 p 

large — + -j i ~ sin 3 , and the SDP (±ir) and L-L* form a closed contour 
which allows us to use the Cauchy Residue Theorem. The pole singu- 
larities contribution yields 


-jk (s+s ' ) 


(A-58) 


tP = i- 


1_ e 

4tt 


1 - 


2ss'sin 3, 


y- (l+cos£ ) 

(s+s') 2 p 


COSNir X 


(s+s' ) 2 -2ss'sin 2 3 0 (l+cos5 p ) 


X H(ir- | -S+2nNir | ) , 

where 

S p = - 3 + 2nN , N = 0, 1, 2, 
and |5 | < tt . 

By making the following approximation 


(A-59) 


2ss'sin 3. 

-jk (s+s ' ) 1 *r- (1+cosc) 

(s+s 1 r 


f ss'sin 2 3. 

-jk (s+s ' )-{l *2. (1+cosC) 

I ( S+S ' ) 


% e 


t . 


the saddle point contributions to Eq. (A-57) can be obtained in a 
manner similar to that for cylindrical wave illumination. Thus 


(A-60) 


'U 



S 1 

s(s+s 1 ) 


e 


-jks 


where d^ is given in Eq. (A-24) with the exception that k = 
in this case. 


kss ‘ 
s+s 1 


sin 2 3. 
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The integral I (p”,z ;p~* ,z' ;e) is evaluated in the same way by 
Pathak and Kouyoumjian[20]. Combining all the results we got thus 
far, and recognizing that the pole contribution for the case N = 0, 
e = <(>-$' yields the direct incident field, we can determine that 

C = 1. 

Thus the diffracted field component of U^s”,?' ) can be written 


(A-61 ) Ujj(s.s') * u 1 Da | s |^ s , - y e~ jks 


where D a are given in Eqs. (A-37a) and (A-37b) with the large parameter 


kss Y _ • 2 0 . 

K " ilV Sln V and 


U 1 * e 


-jks' 


4ttS 1 


86 



REFERENCES 


1. Chu, T.S., "Surface Wave Diffraction and Its Relationship to 
Surface Wave Antennas," Ph.D. Dissertation, The Ohio State 
University ElectroScience Laboratory, Department of Electrical 
Engineering, 1960. 

2. Chu, T.S., Kouyoumjian, R.G., Kara! , F.C., and Karp, S.N., 

"The Diffraction of Surface Waves by a Terminated Structure in 
the Form of a Right-angled Bend," IEEE Trans, on Antenna and 
Propagation, Vol . AP-10, No. 6, November 1962. 

3. Felson, L.B., "Field Solution for a Class of Corrugated Wedge 
and Cone Surfaces," Polytechnic Institute of Brooklyn, Microwave 
Res. Inst., Tech. Memo., No. 32, July 1957. 

4. Harrington, R.F., Time Harmonic Electromagnetic Fields , New 
York, McGraw-Hill, 1961, Chapters V and VI. 

5. Hutchins, D.L. and Kouyoumjian, R.G., "Asymptotic Series 

Describing the Diffraction of a Plane Wave by a Wedge," Report 
2183-3, 15 December 1969, The Ohio State University Electro- 
Science Laboratory, Department of Electrical Engineering; 
prepared under Contract AF 19(628)-5929 for Air Force Systems 
Command. (AD 699 228)(AFCRL-69-0412) . (See also: Hutchins, 

D.L. and Kouyoumjian, R.G., "A New Asymptotic Solution to the 
Diffraction by a Wedge," URSI Spring Meeting, Ottawa, Canada, 
1967.) 

6. Hwang, Y.M., Kouyoumjian, R.G. and Pathak, P., "The Radiation 
from Slots in Truncated, Dielectric-covered Surfaces," 1971 
G-AP International Symposium, Los Angeles, California, 

September 1971. 

7. Karal , F.C. and Karp, S.N., "Diffraction of a Plane Wave by 
a Right Angled Wedge Which Sustains Surface Waves on One 
Surface," New York University, Inst. Math. Sci . , Div. Electro- 
magnetic Res., Res. Rep. EM-123, 1959. 

8. Karal, F.C., Karp, S.N., Chu, T.S. and Kouyoumjian, R.G., 
"Scattering of a Surface Wave by a Discontinuity in the Surface 
Reactance on a Right-angled Wedge," Commun. Pure Appl . Math., 
Vol. 14, February 1961, pp. 35-48. 


87 



9. Karp, S.N. and Karal , F.C., "Vertex Excited Surface Waves on 
One Face of a Right Angle Wedge," New York Univ., Inst. Math. 

Sci . , Div. Electromagnetic Res., Res. Rep. EM-124, 1959 . 

10. Karp, S.N., "Two Dimensional Green's Function for a Right 
Angled Wedge Under an Impedance Boundary Condition," Commu. on 
Pure and Appl . Math., Vol . XIII, 1960, pp. 203-216. 

11. Karp, S.N. and Karal, F.C., Jr., "Generalized Impedance Boundary 
Conditions with Applications to Surface Wave Structures," EM 

wave Theory Proceedings of a Symp. held at Delft, The Netherlands, 
September 1965. 

12. Keller, J.B., "Geometrical Theory of Diffraction," Journal of 
the Optical Society of America, Vol. 52, February 1962, pp. 116- 
130. 

13. Kouyoumjian, R.G., "Asymptotic High-Frequency Methods," 

Proceedings of IEEE, Vol. 53, No. 8, August 1965, pp. 864-876. 

14. Kouyoumjian, R.G. and Pathak, P., "The Dyadic Diffraction 
Coefficient for a Perfectly Conducting Edge Structure," 
presented at URSI Fall Meeting, Los Angeles, California, 1971. 

15. Lewy, H., "Waves on Sloping Beaches," Bull. Amer. Math. Soc., 

Vol. 52, 1946, pp. 737-773. 

16. MacDonald, H.M., Electric Waves , Cambridge University, Press, 
Cambridge, 1902, pp. 186-198. 

17. Malyughinets , G.D., Dokl . Akad. Nauk SSSR, 121, No. 3, 1958, 
p. 436. 

18. Malyughinets, G.D., Annalen der Physik, 6, No. 1-2, 1960, 
pp. 107-112. 

19. Mohsen, A. and Hamid, M.A.K., "Diffraction by a Dielectric-loaded 
Wedge," Radio Science, Vol. 8, No. 1, January 1973, pp. 71-80. 

20. Pathak, P. and Kouyoumjian, R.G., "The Dyadic Diffraction 
Coefficient for a Perfectly-conducting Wedge," Report 2183-4, 

5 June 1970, The Ohio State University ElectroScience Laboratory, 
Department of Electrical Engineering; prepared under Contract 
AF 1 9 (628) -5929 for Air Force Systems Command. (AFCRL-69-0546) 


88 



21. Peters, A.S., "Water Waves over Sloping Beaches and the Solution 

of a Mixed Boundary Valve Problem for (a 2 _ = 0 in a 

Sector," Commu. on Pure and Appl . Math., 5, 1952, pp. 87-108. 

22. Senior,- T .B. A. , "Diffraction by an Imperfectly Conducting Wedge," 
Commu. on Pure and Appl. Math., 12, 1959, pp. 337-372. 

23. Senior, T.B.A., "Impedance Boundary Conditions for Statistically 
Rough Surfaces," Appl. Sci . Res., Section B, Vol . 8, 1961, 

pp. 437-462. 

24. Stoker, J.J., "Surface Waves in Water of Variable Depth," 

Quart. Appl. Math., Vol. 5, 1947, pp. 1-54. 

<t i - . 


ij i 


,y U.S. GOVERNMENT PRINTING OFFICE- 1974— 739-160 ,'122 


< l : .'ft 


89 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE *300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute . . . jo the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS. Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS- 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 





